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Abstract 

Aims. We have investigated the differences in apparent opening angles between the parsec-scale jets of the active galactic nuclei 
(AGN) detected by the Fermi Large Area Telescope (LAT) during its first three months of operations and those of non-LAT-detected 
AGN. 

Methods. We used 15.4 GHz VLB A observations of sources from the 2 cm VLB A MOJAVE program, a subset of which comprise 
the statistically complete flux density limited MOJAVE sample. We determined the apparent opening angles by analyzing transverse 
jet profiles from the data in the image plane and by applying a model fitting technique to the data in the (;(, v) plane. Both methods 
provided comparable opening angle estimates. 

Results. The apparent opening angles of y-ray bright blazars are preferentially larger than those of y-ray weak sources. At the same 
time, we have found the two groups to have similar intrinsic opening angle distributions, based on a smaller subset of sources. This 
suggests that the jets in y-ray bright AGN are oriented at preferentially smaller angles to the line of sight resulting in a stronger 
relativistic beaming. The intrinsic jet opening angle and bulk flow Lorentz factor are found to be inversely proportional, as predicted 
by standard models of compact relativistic jets. If a gas dynamical jet acceleration model is assumed, the ratio of the initial pressure of 
the plasma in the core region Pq to the external pressure Pext lies within the range 1 . 1 to 34.6, with a best fit estimate of Pq /Ptxt ~ 2. 
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1. Introduction y-ray sources, 116 of which were associated with high galactic 

latitude i\b\ > 10°) AGN (the LAT Bright AGN Sample together 

The EGRET telescope onboard the Compton Gamma Ray wi th 10 lower-confide nce associations: LBAS-ext), as discussed 

Observatory dHartman et alj 1 19991) provided detections of y- bv lAbdo et al. l ('2009a'). Their VLBI radio source identifications 

ray emission from many extragalactic poi nt sources, rt i ost of were reported by Kovalev (2009). The LAT detections of flaring 

which were identified with blazars (e.g., iMattox et"an 120011: AGN have also triggered a nu mber of multiwave length (from 

ISowards-Emmerd et al.ll2003[) . The latter term generally refers radio to y-ray) campaigns (e.g.. lAbdo et alj|2009bl) . 

to objects classified as flat-spectrum radio-loud quasars and BL The Fermi era has akeady heralded a number of impor- 

Lacs. Although blazars comprise only a few per cent of the over- tant results that link the y-ray emission and radio properties 

, an AGN population, they dominate the extragalactic high-energy of AGN. It has been shown that LAT-detected blazars are 

• sky. A systematic comparison of parsec-scale radio jet structure brighte r and more luminou s in the radio domain at parsec 

with y-ray emission in AGN has hitherto been problematic and scales ( Kovalev et alj l2559h , have higher apparent jet speeds 

inconclusive due to the limited sensitivity and angular resolu- (|Lister et al. 2009g). The LAT-detected blazars also seem to have 

tion of the EGRET instrument. For example, a recent analysis higher Do ppler factors determine d from die mm-wavelength 

of a large 6 cm VL BA Imaging and Polarimetry Survey (YIPS; variability (ISavolainen et al.ll2009h . The y-ray and radio flares 

[Taylor et al.| |2007|) has hinted fliat EGRET-detected blazars are found to appear within a typical timescale of up to a few 

might have larger than average jet opening angles, but very poor months and are likely associated with the parsec-scale radio core 

statistics (only four EGRET/VIPS sources) did not allow the au- (Kovalev et aL 2009)- In this Letter we continue investigating the 

thors to draw a firm conclusion. relations between y-ray brightness and the properties of parsec- 

On 2008 June 11, the Fermi Gamma-Ray Space Telescope scale radio jets from the LBAS-ext list as probed by interfero- 

(previously known as GLAST) was successfully launched by metric observations with the VLBA. Hereafter we refer to the 

NASA with the Large Area Telescope (LAT), a successor to Fermi LAT 3-month bright source list dAbdo et alj [ 2009al) when 

EGRET, onboard. The LAT began operating in August 2008 and using the term LAT-detected. 

has provided y-ray observations with greatly improved sensitiv- Throughout this letter, we use the term "core" as the appar- 

ity, superior angular resolution, large field of view , and wide en- ent origin of AGN jets that commonly appears as the bright- 

ergy range from about 20 MeV to over 300 GeV ( lAtwood et alj est feature in V LBI images of blazars (e.g., Lobanov 1998|; 

|2009|) . During the first three months of science operations, the lMarsche3l2008[) . We use the ACDM cosmological model with 

LAT has made 205 high-confidence (> lOcr) detections of bright Hq -71 km s"' Mpc"', f2„, = 0.27, and Q.^ = 0.73. 
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Figure 1. Left: Naturally weighted total intensity CLEAN image of 1641+399 at 15 GHz observed on 2009 February 25. The jet 
axis is approximated by two straight lines in position angles of -90° and -50°. The contours are plotted at increasing powers of 2, 
starting from 0.4 mJy beam"'. The peak flux density reaches 4770 mJy beam"'. The FWHM of the restoring beam is shown as a 
cross in the lower left corner. Middle: Deconvolved FWHM transverse size of the jet along its axis. Right: Apparent opening angle 
of the jet along its axis. 



2. The radio data and source samples 

The MOJAVE program dLister et alj l2009a[) is a long-term 
VLBA program to study the structure and evolution of ex- 
tragalactic relativistic radio jets in the northern sky. The full 
monitoring list currently consists of more than 250 sources, 
and includes a statistically complete, flux-density limited sam- 
ple of 135 AGN (hereafter referred to as MOJAVE-1). All the 
MOJAVE- 1 sources have J2000 declination 6 > -20° and a 
15 GHz VLBA correlated flux density Sc > 1.5 Jy (2 Jy for 
6 < 0°) at any epoch between 1994.0 and 2004.0. The weaker 
radio blazars [Sc > 0.2 Jy) detected by Fermi extend the com- 
plete MOJAVE-1 sample to MOJAVE-2. 

Although we have analyzed a large portion of the MOJAVE- 
2 sample, we focus our statistical analysis on a MOJAVE-1 
blazar subset (101 quasars and 22 BL Lacs), which is a dominant 
fraction of the sample (91%). We excluded eight radio galax- 
ies and four optically unidentified objects from the MOJAVE-1 
sample. Taking into account the fact that the LBAS-ext list is 
restricted in galactic latitude {\b\ > 10°), due to the decrease of 
the LA T detection sensitivity at low galactic latitude (lAbdo et al.l 
I2009ah . seven quasars and one BL Lac were excluded from the 
analysis as having \b\ < 10°. Of the remaining 115 MOJAVE-1 
blazars, 29 are in the LBAS-ext list, corresponding to an overall 
y-ray detection rate of 25%. It is worth noting that the detection 
rate for BL Lacs is substantially higher (48%) compared to that 
for quasars (20%) for the MOJAVE-1 sample. 

3. Results 

The opening angle of a jet in the sky plane can be estimated 
either by analyzing transverse jet profiles in the image plane, or 
by model fitting the interferometric visibility data in the (u,v) 
plane. To increase the robustness of our conclusions we used 
both methods and compared the results. 

3.1. Opening angles from the image plane 

We used the 15 GHz naturally weighted MOJAV^] VLBA im- 
ages from the most recent epoch available. This covers the pe- 

* http://www.physics.purclue.edu/MOJAVE 



riod 2008-2009 for all but two sources. The opening angle of the 
jet was calculated as the median value of 0- = 2 arctan[0.5((i~ - 
b^y^^/r], where d is the full width at half maximum (FWHM) of 
a Gaussian fitted to the transverse jet brightness profile, r is the 
distance to the core along the jet axis, b^p is the beam size along 
the position angle (p of the jet-cut, and the quantity {d^ - b^y^^ 
is the deconvolved FWHM transverse size of the jet. The direc- 
tion of the jet axis was determined using the median position 
angles of all jet components over all the epochs from model fit- 
ting (see next section). The slices were taken at 0. 1 mas intervals 
starting from the position of the VLBI core and continuing up to 
the region in which the jet either substantially curved or became 
undetectable. The ridge Unes for 15 MOJAVE-1 and 6 MOJAVE- 
2 sources with notably bending jets were approximated by two 
straight lines. Opening angle values were calculated using only 
those slices that had a peak of the fitted Gaussian larger than 4 
times the rms noise level of the image. In Fig.[T]the 15 GHz total 
intensity map of 3C 345 (1641H-399), together with the decon- 
volved size and opening angle of the jet as a function of angular 
distance to the core is shown as an example. 

The distributions of measured opening angle are shown in 
Fig. |2] The distribution for LAT-detected sources is narrower 
and missing the small opening angle sources if compared to 
non-LAT-detected sample. A Kolmogorov-Smirnov (K-S) test 
indicates a probability of only p - 0.019 for these two sam- 
ples being drawn from the same parent population. If we add 27 
additional LAT-detected blazars from the extended MOJAVE-2 
sample the confidence level increases to 99.9% {p < 0.001) and 
the mean value to 22^9 + 1^5 (Fig.|2l bottom panel). The jet-slice 
opening angles are on average larger for the 29 LAT-detected 
MOJAVE-1 blazars in comparison with the 86 non-LAT-detected 
ones, with a mean value of 21?4 ± V.S versus 18?0 + l^l, re- 
spectively. A Student's T-test confirms that the means are dif- 
ferent at 96.1% confidence level for the complete MOJAVE-1 
sample and at 99.4% confidence level for the MOJAVE- 1,2 sam- 
ple. The relatively small diff'erence in mean values is due to a 
significant number of y-ray weak sources with very large ap- 
parent opening angles. In fact, this is compatible with the find- 
ings of [Savolainen et al. (2009), who show that there should be 
a number of y-ray weak sources that are viewed nearly end-on 
and which therefore have very large apparent opening angles. 
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Our results support this. We note that the TANAMI group, on 
the basis of VLBI analysis of a Southern hemisphere sample of 
Fermi detected AGN jets, has also found indications of prefer- 
entially wider apparent opening angles in LAT-detected sources 
(R. Ojha et al., in preparation). The calculated values of the ap- 
parent opening angles are listed in Table iQ. 

We have studied the effects of a poss ible bias by BL Lac ob- 
jects, as discussed bv lAbdo et al. I (l2009ah . on our results by com- 
paring the apparent opening angle distributions of 22 BL Lacs 
and 101 quasars from the complete MOJAVE-1 sample. The K- 
S test indicated no significant difference (p - 0.740). Another 
check was performed by excluding the BL Lacs from our anal- 
ysis. In this case, the small set of LAT-detected MOJAVE-1 
quasars (A^ - 19) did not produce a statistically significant re- 
sult. However, comparison of the apparent opening angle distri- 
bution for LAT-detected versus non-LAT detected quasars in the 
larger MOJAVE-l+MOJAVE-2 sample led to the following. The 
distributions are significantly different according to both the K- 
S and the Student's T-tests (pk-s = 0.006, /^j-test = 0.017). We 
conclude therefore that the presence of BL Lacs does not bias 
our analysis. 

At the same time, we found no significant correlation be- 
tween the apparent jet opening angle and average F ermi LAT 
100 MeV - 1 GeV photon flux dAbdo et al.ll2009al) . since the 
latter depends on the distance to a source. However, we would 
expect to see a positive correlation between the apparent open- 
ing angle and y-ray luminosity, when the y-ray energy fluxes 
become available in the next Fermi data release. 



3.2. Opening angles from the {u, u) plane 

The estimates of the opening angle can also be derived from the 
data in the («, u) plane using a model fitting approach. The ob- 
served brightness distribution of each source was modeled by a 
limited number of two-dimens ional Gaussian c omponents using 
the "modelfit" task in Difmap jShepherdlll997l). Th e parameters 
of the models are tabulated bv lLister et al.l (12009 bh . We used all 
available 2336 model components for 115 MOJAVE-1 blazars 
with galactic latitude \b\ > 10°. After measuring the opening an- 
gle for each jet component, we averaged them to get one value 
for each epoch, and then averaged over the epochs to obtain a fi- 
nal estimate. The statistical results of the model-fit opening angle 
distribution analysis are in good agreement with those from the 
transverse jet profile method, confirming that our conclusions 
are robust. The median value of the ratio of apparent opening 
angle derived from jet-cut analysis to that from the model fitting 
method is 0.9. The apparent opening angle distributions of LAT- 
detected and non-LAT-detected MOJAVE-1 blazars are different 
at 99% confidence level, according to a K-S test. 



3.3. Intrinsic opening angles 

The on-average wider apparent opening angles in LAT-detected 
sources can be a consequence of (i) their smaller viewing 
angles (due to a projection effect) and/or (ii) their intrinsi- 
cally different properties, for example, the presence of spine- 
sheath stru ctures like those obse rved in the y-ray bright quasar 
1055-H018 (lAttridgeetal.lll999[). TeV BL L ac object MrkSOl 
dPushkarev et al.' 1 20051: iGiroletti etani2008l) and radio galaxy 
M87 (Kovalev 2003)- 

To distinguish between these two possibilities, we have de- 
rived the values of the viewing angle and the bulk Lorentz fac- 
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Figure 2. Distributions of the apparent opening angle from jet- 
cut analysis for 86 non-LAT-detected {top panel, unshaded), 29 
LAT-detected (middle and bottom panel, dark gray) MOJAVE-1 
blazars, and 27 LAT-detected {bottom panel, light gray) addi- 
tional MOJAVE-2 blazars. 



tor r for MOJAVE-1 blazars (see Table 1) using the following 
relations; 



9 - arctan ■ 



app "^var 



r = 



i2 

app 



25v; 



Table 1 is available in electronic form at |http://www.aanda.org| 



where jS^pp is the fastest measured radial, non-acceler ating appar- 
ent jet speed from the MOJAVE kinematic analysis dLister et all 
'2009bl) and (5var is the variability Doppler factor from the 
Metsahovi AGN monitoring program (Hovattaet al. 2009)- The 
overlap of the MOJAVE and Metsahovi programs comprises 56 
blazars with measured speeds and Doppler factors. Indeed, the 
viewing angles of 21 LAT-detected sources turned out to be 
sUghdy smaller, with mean value of 3?6 + 0?4 vs. 5°.l ± 1^3 for the 
non-LAT-detected sources. However, the difference is not statis- 
tically significant for this small sample (56 sources compared to 
115 in the apparent opening angle analysis) containing a mix- 
ture of quasars and BL Lacs. If only quasars are considered, 
a marginally significant difference {p - 0.06) in viewing an- 
gle distributions i s found between LAT-d etected and non-LAT- 
detected sources dSavolainen et al.|[2009l) . In addition, we have 
found an indication for BL Lacs to have on-average wider intrin- 
sic opening angles (2°4 ± 0°6) than those of quasars (1°2 + 0^1). 
The corresponding distributions are different at confidence lev- 
els of 94.6% according to the K-S test, the average values differ 
with a 96.1% confidence according to the Student's T-test. 

To test the second scenario, the intrinsic opening angles 
Q'int = Capp sin© were calculated for the 56 sources (Table 1). 
The estimated values of this parameter range up to 8° (Fig.|3]l. A 
K-S test indicated no significant difference {p - 0.797) between 
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Figures. Intrinsic opening angle vs. Lorentz factor for 56 jets. 
Two sources with F, Cint values of (45,0.96) and (65,0.53) re- 
spectively are beyond the plot limits. The solid line shows the 
median curve fit with the assumed relation a - p/F, where p 
is a constant. Filled circles correspond to LAT-detected sources, 
while open ones correspond to non-LAT-detected sources. The 
dotted curves represent relationships between the opening angle 
and Lorentz factor as predicted by the gas dynamical model for 
different values of the parameter ^. 



the samples of LAT-detected and non-LAT-detected sources, 
suggesting that the established systematic difference in apparent 
opening angles is most probably the result of projection effects, 
i.e., the y-ray bright jets are aligned closer to our line of sight. 

We have also analyzed the observed dependence between 
the intrinsic opening angle and Lorentz factor (Fig. [3]), which 
are expected to be inversely proportional according to simple 
hydro dynamical models of relativistic jets (Blandford & KoniEj] 
Il979l) . Both the gas dynamical m odel ODalv & Marscher., 198^5 ^ 
and magnetic acceleration models dKomissarov et al.ll2007h also 
predict this relation. The observed dependence was fitted assum- 
ing a relation aim = p/F with the coefficient p left as a free 
parameter. The best fit value of p was found to be 0.26 rad by fit- 
ting the median curve. (Fig.[3j solid line). Even though the flow 
Lorentz factor and the viewing angle (which is needed to deter- 
mine ffint), are both calculated from the same observables (ySapp 
and 6), a Monte Carlo simulation shows that the apparent nega- 
tive correlation between aim and F is not due to this degeneracy, 
but is a genuine effect. 

In the gas dynamical model of compact relativistic jets sug- 
gested bv iDalv & Marsche^ ( Il988h . the opening angle of a jet 
is a function of the Lorentz factor and a ratio of the external 
pressure Pgxt to the initial pressure Pq of the plasma in the core 
region, ^ - V^ext/A)- We applied this model for different values 
of ^ (Fig.|3] dotted curves) and were able to constrain the param- 
eter ^ to lie within a range [0.17,0.96] with the best fit estimate 
f = 0.67, corresponding to a range of [ L 1 , 34.6] for Pq / Pext with 
the best fit estim ate Pp/Pext ~ 2. Ou r results confirm those ob- 
tained earlier by iJorstad et al.l ( l2005l) using a smaller sample of 
15 blazar jets. 



4. Summary 

We have measured the projected jet opening angles on par- 
sec scales for 115 blazars (29 LAT-detected and 86 non-LAT- 
detected) from the complete flux-density limited MOJAVE-1 
sample and for 27 additional LAT-detected sources monitored 
by the MOJAVE program. The apparent opening angles for y- 
ray bright sources are on average larger than those in y-ray weak 
ones, while the intrinsic opening angle distributions based on 
smaller samples are statistically indistinguishable. We interpret 
this as an evidence for y-ray bright blazars to have preferentially 
smaller viewing angles and, consequently, stronger relativistic 
beaming, which boosts emission in both, the y-ray and radio 
bands. This conclusion is consistent with recently obtained re- 
sults that show LAT-detect ed blazars to be bri ghter and more 
lumi nous in radio doma in dKovalev et al.l 120091) . to have faster 
j ets (Lister et al."2009c') and higher variability Doppler factors 
dSavolainen et al. 2009). There is an indication for BL Lac ob- 
jects to have on-average wider intrinsic opening angles than 
those of quasars. 

The intrinsic opening angle and Lorentz factor are found to 
be inversely proportional for a sample of 56 AGN jets in agree- 
ment with theoretical predictions of simple gas dynamical and 
magnetic acceleration models. The best approximation of the in- 
ferred relation is achieved assuming a ratio of internal and exter- 
nal jet pressure Po/Pext ~ 2 in gas dynamical models of compact 
relativistic jets. 
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